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and subsequently apply these operations at finite temperatures to investigate the effect 
of decoherence on operational fidelity. 
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1. Introduction 

Quantum computation is an extremely active area of research, with potential realizations 
ranging from solid-state to atomic and optical systems Pfl |2j. Following seminal 
work on ideal quantum computation, the focus of research has shifted in recent 
years to practical issues of implementing scalable quantum gates in real systems in 
the presence of noise, decoherence, and imperfect controls. One significant issue 
is the fundamental process of decoherence caused by coupling to the environment. 
Another obstacle is noise that, in any practical situation, exists in external control 
fields which are normally required to perform the logical operations of quantum 
information processing 3 . In order to ensure fault-tolerant operation of quantum 
registers, errors, including those associated with environment ally- induced decoherence 
and control field noise, must be minimized below a specified threshold [IJ. Since 
external controls are essential for logical operations, it is reasonable to assume that 
practical implementations of quantum computation, especially in relatively complex 
physical systems, can be greatly facilitated by the powerful methods of optimal control 
[Slinj. Closed-loop adaptive laboratory control is of special importance because of its 
inherent robustness to noise and suitability for managing complex physical systems 
El HI El EE IHIl EH EE2j • Furthermore, intense ultrafast control fields, which are essential 
for accelerating quantum operations, thereby diminishing the effect of decoherence, are 
especially suited for closed-loop laboratory control. Applying the methods of optimal 
control to various realizations of quantum computation was recently considered in a 
number of works HS1 HH [13 HH1 HZl HH1 HH1 123 |2H 1221 

The majority of schemes proposed for physical realizations of quantum computation 
assign a single physical level to represent each of the logical basis states (|0) and |1)) of 
a qubit. In such a situation, each logical operation is represented by a unique physical 
transformation pQ. In this paper we consider a more general situation, where each logical 
basis state is encoded into a subspace of multiple physical levels. This generalization 
leads to the logical equivalence of different physical states, i.e., different multilevel 
superpositions can represent the same logical state of a qubit. Correspondingly, with 
multilevel encoding (MLE) there will exist different physical transformations which 
will also be logically equivalent, i.e., these transformations represent the same logical 
operation. 

An important question is whether the use of MLE yields any practical advantages 
in comparison to the usual single-level encoding (SLE). The following analysis indicates 
that MLE is a natural choice when the need to accelerate quantum operations requires 
the bandwidth of control fields to exceed the separation between energy levels of a 
quantum register. For systems with a dense energy spectrum, the factor determining 
individual level addressability is the ratio between the bandwidth, AQ, of the external 
control field and the characteristic energy-level spacing, Au, of the physical system that 
realizes the register. In the limit of slow control, i.e., AQ < Auj, logical basis states 
can be encoded into single physical levels and MLE is not necessary. While such slow 
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controls may be satisfactory for proof-of-principle experiments, practical realizations of 
quantum registers in many cases will require much faster operations and correspondingly, 
faster controls, with AQ 3> Au. In such a case, one must take into account the 
various transitions between the multiple states driven by a control field. Encoding 
each logical basis state into a subspace of multiple physical levels circumvents this issue, 
since addressability of individual levels within each subspace is not necessary due to the 
principle of logical equivalence developed in this work. 

It is clear from the arguments above that the importance of MLE depends on the 
specific physical realization of a quantum register, so it will be helpful to consider specific 
examples of systems which may be suitable for realizations of MLE. One such system 
is a molecule in which the logical basis of a qubit is realized by vibrational states on 
the ground and/or excited electronic surfaces |23| I24j . For example, typical vibrational 
periods in homonuclear alkaline dimers range from 300 to 500 fs [25 . If the system 
is driven by ultrafast laser pulses of durations ranging from 20 to 100 fs (which are 
currently successfully used in closed-loop adaptive control experiments with molecular 
systems [211 123 I2E1 I2H ELU EH E2]), many transitions between vibrational levels will be 
driven simultaneously. In such a case, it can be useful to encode each logical basis state 
into a set of multiple molecular levels. 

Another interesting example is a system whose internal levels are used for the 
logical basis while it is spatially confined in an external potential (e.g., an ion in an 
electromagnetic trap or an atom in an optical lattice). In the majority of schemes 
proposed for quantum computation with trapped ions, the gate time is limited by 
the trap frequency, mainly due to the need to spectroscopically resolve the motional 
sidebands jH3 OH IHH HSU EZ|- This is a serious limitation because typical ion-trap 
frequencies are relatively low, ranging from 100 kHz to 100 MHz |38| l3*Uj . External 
controls which are faster than the vibrational frequency of a trapped qubit will inevitably 
excite multiple motional levels. This effect will be especially important outside of the 
Lamb-Dicke limit, i.e., when multiple sidebands are excited by a laser field interacting 
with a trapped ion. Encoding each logical basis state into multiple motional levels of a 
trapped particle can be a useful approach in situations when fast controls are desired. 

An idea related to the concept of MLE was recently proposed by Garcia- Ripoll et al. 
[4*0] , who considered a two-qubit operation acting on a pair of trapped ions. The essence 
of this proposal is to apply coherent laser controls which leave the initial motional state 
unchanged at the end of the sequence of control pulses, thereby allowing for operations 
which are much faster than the trap frequency. MLE allows for more general operations 
(not just the identity operation) on the motional levels which constitute the encoding 
subspaces for a trapped particle. However, MLE may not be directly applicable to 
popular schemes where the common motional mode of a number of trapped ions is 
additionally employed to entangle the internal atomic states of different qubits. At 
this point it is not yet clear which systems will be most suitable for MLE in practical 
realizations of quantum computation, but the importance of MLE will increase in line 
with the acceleration of quantum operations. 
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Searching for an optimal control field to perform a logical operation on a multilevel 
system can benefit from MLE due to the logical equivalence of different physical 
transformations. Also, the multiplicity of logically equivalent physical transformations 
can facilitate finding optimal operations with improved robustness against noise in 
external controls [H]. Moreover, multiple physical transformations which realize the 
same logical operation may be used interchangeably throughout the computation, 
providing greater flexibility in the laboratory implementation of controls. 

Perhaps the most important property of MLE is its utility as a practical instrument 
to reduce the effects of mixing in the initial quantum state (e.g., due to thermal 
excitations) and decoherence throughout the computation. In many practical situations, 
the initialization of a quantum register requires cooling a system to the ground state 
since the impurity of the initial state due to thermal excitations reduces quantum- 
computational fidelity. This cooling is often a difficult and slow process, therefore a 
practical alternative is desirable. Another significant source of errors is decoherence 
during logical operations. A possible way to overcome these problems is, once again, 
by encoding the logical basis into subspaces of multiple physical levels. Due to the 
principle of logical equivalence, mixing and decoherence within each encoding subspace 
do not affect coherence between the two logical basis states and therefore do not reduce 
computational fidelity. 

This paper is organized as follows. In section El the concept and mathematical 
formalism of multilevel encoding are developed for a single qubit. Logical equivalence 
resulting from MLE is defined for physical states and transformations, and equivalence 
classes of logical states and operations are introduced. In section EJ the MLE 
formalism and corresponding logical equivalence are used to construct single-qubit 
unitary operations and derive their general form (a tensor-product structure). Section 0] 
extends the MLE formalism to mixed states and non-unitary operations. In particular, 
this allows for logical equivalence between pure and mixed states, and between unitary 
and general non-unitary (permitting decoherence) quantum operations, provided that 
full coherence exists between the encoding subspaces. In sectional we demonstrate the 
formalism and some of its advantages by performing a numerical simulation of closed- 
loop optimal control of single-qubit operations for a model system based on the electronic 
and vibrational levels of an alkaline dimer, and applying the resulting operations at finite 
temperatures. In sections HI and [71 the MLE formalism is generalized for logical states 
and operations of multiple qubits. As in the single-qubit general tensor-product 

form for multi-qubit unitary operations is obtained. Section |H] concludes the paper with 
a brief summary of the results, open questions, and future directions. 

2. Concept and Formalism 

For a single qubit, MLE involves partitioning the energy levels of a given quantum 
system into two distinct subspaces of equal dimension, the "encoding subspaces" S 
and Si. Let \a) denote any element of So such that (a\a) = 1 and |6) denote any 
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element of Si such that (6|6) = 1, i.e., 

n 

i°) = ^2 a i\xi) = ( 



T 

a± ■ ■ ■ a n 



i=i 



such that (a\a) = |<2j| 2 = 1 (1) 

8=1 

and 

\b) = J2bM = (h •••fen) 

i=l 

n 

such that (b|b) = ^|6i| 2 = 1. (2) 

i=i 

Here n is the dimension of the encoding subspaces and and {|0j)} are orthonormal 

bases that span S and <Si, respectively. Correspondingly, (xi\Xj) = = $ij an d 

(Xil^j) — 0. It is important to emphasize that for each subspace (So or <Si) there exists 
an infinite set of states that satisfy the respective criteria of (J2J) or (J2J). Within this 
framework each logical basis state is encoded into one of these subspaces (figure HJ). 
Therefore, the Hilbert space of an encoded qubit (which is two-dimensional in the SLE 
representation) corresponds to an expanded (and higher-dimensional) Hilbert space of 
physical levels. In this expanded Hilbert space the logical basis states are expressed as 

|0) = |0) L ® I a) = ••■ a n Q_^_0,) T (3a) 

n 

|1> = |1) L <8> 1 6) = (^^O h ■■■ b^, (3b) 



where |0)l = (1>0) T and = (0, 1) T are the "logical components" of the logical basis 
states and |o) and |6) are the "encoded components." It follows that (0|0) = (1|1) = 1 
and (0|1) = 0. Thus, any superposition state \ip) can be written as 

|V) = Cb|0) + Cl |l) 

= c |0) L ® \a) +d|l) L ® \b), (4) 

where is a vector of dimension 2n and |co| 2 + |ci| 2 = 1 so that (ip\ip) = 1. Throughout 
this work we employ the following abbreviated notation: 

(5) 

v h> ) 

where a = c \a) and (3 = ci [ 6) . 

There are many different physical situations in which a tensor-product structure 
(TPS) in the form of (jlj) can be realized [12]. An important point is that in the case 
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Figure 1. A sample scheme that illustrates the difference between single-level and 
multilevel encoding of logical basis states of a qubit. 



of MLE this TPS is virtual since it does not arise from spatially separated degrees 
of freedom, instead it is an abstract partitioning of the Hilbert space. This essential 
feature makes MLE physically quite different from the mathematically similar structure 
of quantum error correction (QEC) In the case of QEC one deals with a number 
of distinct physical systems (a qubit and ancillae) which become entangled via a set of 
controlled interactions between them. However, in the context of MLE we deal with a 
single quantum system for each qubit (even if this system is in fact composed of many 
particles, e.g., a molecule), and use of the tensor product is a convenient mathematical 
description. The Hilbert space of a molecule or atom is just one example of a state space 
in which a virtual TPS can be built. Other, very general mechanisms of producing a TPS 
come from superselection rules and the representation theory of algebras of observables 
|42j . Therefore, there is a plethora of possibilities for building an MLE structure: from 
single atoms to molecules, to complex multiparticle systems. 

As seen from (jHJ) and (JH), different physical states, i.e., different choices of \a) and 
1 6), can represent the same logical state. In order to rigorously formulate this principle 
of logical equivalence, we introduce a set of "equivalence assay" (EA) operators: 

1 n 

A* = 2^(1^1 + 1^1)' ( 6a ) 

i=l 
1 - 

i=l 
1 - 

a, = gEflx*) foi-ifcXfci)- ( 6c ) 
i=i 

The EA operators can also be rewritten in the matrix form 

A r = A r <g> l n = -oy ® t n , r = {x,y,z}, (7) 

where {a x , a y , cr z } are the three Pauli matrices, {A^, A^, X z } are the 2x2 matrices 
(A r = \a r ) which represent the three generators of the SU(2) group, and l n is the 
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n x n identity matrix. Expectation values of the EA operators of an arbitrary state 
are 

MA*^) = Re(c* Cl (a|6)) = Re (a f /3) , (8a) 
(VIA^IV) = Im(c* Cl (a|&)) = Im , (8b) 

(^|A^> = ^(|c | 2 (a|a) - | Cl | 2 (fe|6>) = ~ (|c | 2 - | Cl | 2 ) . (8c) 

Two physical states, \ip) and \ip'), are logically equivalent, which is denoted as ~ |?//), 
if they satisfy 

ty\KW) = W\&rW), r = {x,y,z}. (9) 

In addition, the states must be normalized ((4>\4>) — (^\^-2 <8> J-nl^) — 1 f° r a U l^))- 
Logical equivalence resulting from MLE is analogous to equality in the context of SLE. 
Since the encoded component of all EA operators is the identity operator l n , logical 
equivalence, as defined by (J5)l. allows for significant freedom in the choice of the encoded 
components of a state without making it logically different. On the other hand, the 
logical components of the EA operators, oy, together with the identity operator 1 2 , 
span the space of all 2 x 2 matrices, ensuring that all equivalent states are logically the 
same, i.e, that no logical operation can distinguish between them. 

Logical equivalence is a transitive property, i.e., if \ip) ~ \ip') and ~ W), then 
W) ~ W)- Therefore, all states which are mutually equivalent form an equivalence 
class 8$, which is formally expressed as 

£4 = {w.) : \a) ~ my. (10) 

Correspondingly, every physical state in a given equivalence class represents the same 
logical state. 

Logical equivalence can also be extended to quantum operations. Two operations, 
U and U' are logically equivalent (denoted as U ~ U'), if 

U\ip) ~ U'\il)) for any \tfj). (11) 

Logical equivalence for operations is also a transitive property, i.e., if U ~ U' and 
U ~ U", then U' ~ U" . Therefore, all physical transformations which are mutually 
equivalent form an equivalence class and every transformation in that class represents 
the same logical operation. 

3. Single-Qubit Operations for Multilevel Encoding 

Any 2x2 unitary matrix can be represented as a linear combination of the identity 
matrix II2 and the three Pauli matrices {o~ x ,a y ,a z }. The Pauli matrices (up to a 
numerical factor) are the generators of the Lie group SU(2) and satisfy a set of 
commutation and multiplication properties. Single-qubit logical operations for SLE are 
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represented by elements of the U(2) group. In this section we develop the corresponding 
representation of single-qubit logical operations within the context of MLE. 

Using the partitioning scheme for logical basis states described in section the 
general form of a unitary operation acting on a qubit encoded in the expanded tri- 
dimensional space is 

°-{ A o B n)> 

where each sub-block A, B, C, and D is an n x n matrix representing operations on and 
between the two encoding subspaces, So and S\. Note that any operation U is defined 
up to an overall phase factor. 



3.1. The weak identity operation 

The weak identity operation, [//, is defined by the following property: under the 
conditions of logical equivalence, Uj leaves the equivalence class E^ invariant, i.e., 
UjE^j, = E^, which can be expressed as 

W)HV>> for any 1^). (13) 
Using the general decomposition of (fT2"j). the equation above is recast in the form 



Aa + Bp \ la 
Ca + D{3 J \ p 



(14) 



We find the explicit form of Ui by using the conditions of logical equivalence, (jHJ). The 
condition for A 2 , combined with the normalization property, requires that 

+ /?t£t) (Aa + BP) = a f a, (15a) 

(a t C t + P^D^) (Cat + D/3) = ft (3, (15b) 

which is true if A and D are unitary and B — C — 0. Using these results, we next 
consider the two other conditions for equivalence (A x and A y ): 

j P±p^a, (16) 

which are satisfied if A = D. Therefore, we obtain 

U, = I V ° ° I = II Vo. (17) 



Vr 







where Vq is an arbitrary n x n unitary matrix and denotes the n x n matrix of zeros. 
Therefore, the weak identity Uj is in fact a class of logically equivalent operations: 



Uj = {1 2 ®V : Vq E U(n)}. 



(18) 
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The identity operation l 2n is also a member of this equivalence class (corresponding to 

V = l n ). 

The notion of the weak identity operation clarifies the meaning of logical equivalence 
between different physical states. For any state = (a,j3) T , the equivalence class £</, 
is generated by the set of unitary operations Ui acting on 

= {(UVVqM) :V Q eV(n)}. (19) 

It is impossible to preserve logical equivalence when different transformations act on 
a and j3, e.g., the state \ip) = (Va, V [3) T with V ^ V is logically different from the 
state \ip) = (a,/3) T since the relative phase between |0) and |1) is changed. This is why 
the tensor-product structure of (|T7|) is necessary for maintaining logical equivalence and 
emphasizes why logical equivalence is defined for a superposition state \ip). It would be 
meaningless to define logical equivalence for the individual logical basis states since no 
definite relative phase between |0) and |1) would exist. 

It follows from (fT9~|) that an infinite number of different physical states within an 
equivalence class correspond to the infinite number of different choices for the unitary 
transformation Vq. There also exists an infinite set of equivalence classes themselves. 
If \ipi) = (ai,/?i) T and 1^2) = («2,/?2) T are such that Ui\ip\) 7^ \ip2) for any choice of Vo, 
then the two equivalence classes and S^ 2 are distinct and have no common states. 
Theoretically there is no preference as to what equivalence class to use for MLE of a 
qubit, however, in a practical situation some specific equivalence class will be selected, 
for example, depending on the initial state of the qubit physical system or the technical 
limitations of the available controls. 



3.2. The Pauli matrices 

We will denote the high- dimensional representations of the Pauli matrices {c x ,a y ,a z } 
as {U x ,U y ,U z }. The construction of these operations for a multilevel-encoded qubit is 
similar to that of Ui. In particular, the form of the bit-flip operation U x (the NOT gate) 
under the condition of logical equivalence is 




U *=[ Vl l=^®^i» (2°) 
where V\ is an arbitrary n x n unitary matrix. It is easy to verify that 





As demonstrated in (}2*T]) . U x performs the desired bit-flip operation. Analogously, the 
other two Pauli matrices U y and U z (the latter being the phase-flip operation) have the 
form 

-iV 2 
iV 2 



U y = _ n 2 \=a y ®V 2 (22) 
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u * = ( o 3 -V 3 J = ® ^' ^ 

where V 2 and V 3 are also arbitrary nxn unitary matrices. Under logical equivalence, the 
multiplication rule Oj = iejkiO'k 'i (where {j, k, 1} = {x, y, z} and ejki is the completely 
antisymmetric unit tensor) and the anti-commutation rule {aj,a k } = (j 7^ k) for the 
Pauli matrices are expressed as 

Uj ~ ie jk iU k Ui, {j, k, 1} = {x, y, z) (24a) 
U 3 U k ~-U k U 3 (j^k), (24b) 

respectively. One can easily verify that these relations are satisfied given the form of 
U x , U y and U z and the arbitrary nature of Vi, V 2 , and V 3 . 

3.3. The Hadamard gate 

We denote the high-dimensional representation of the Hadamard gate H as Uh- Since 
Ha z H = a x , Uh must satisfy UhU z Uh\iP) ~ U x \ip) for any \ip). Consider Uh of the 
most general form of (|12jl : 



AV 3 A-BV 3 C AV 3 B-BV 3 D \f a \ I (3 
CV 3 A - DV 3 C CV 3 B - DV 3 D )[ (3 ) ~ 1 a 



(25) 



which is true if A = B = C = —D and \[2A is unitary. Therefore, Uh has the form 

where V4 is an arbitrary nxn unitary matrix. 
3.4- The phase-shift gate 

We denote the high-dimensional representation of the phase-shift gate P(4>) as Up(<fi), 
which has the form 

u P {<P)=[ V * v ° e <+)= p W® v *> ( 27 ) 

where G M and V5 is an arbitrary nxn unitary matrix. 
3.5. Weak commutation 

Under logical equivalence, we can introduce a weak commutation relation. Two 
operators F and G are said to weakly commute if 



FG ~ GF. 



(28) 
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Within the MLE formalism, every logical operation should weakly commute with Uj 
(the weak identity). It is not difficult to see that 17 i weakly commutes with the logical 
operations U x , U y , U z , Uh and Up(4>) described above. Moreover, we show that Ui 
weakly commutes with any unitary operation 17% such that 

C^eU(2)®U(n), (29) 

which is a subgroup of U(2n). In order to prove that 

UjUt ~ E^Z7/, (30) 

first note that any such operation can be written as 17% = £ <8> V, where £ G U(2) and 
V G U(n). Then we find that 

(i)\U\u\KUiU^) = Tr ((£%£) ® (31) 

and 

{^UplKUtU^) = Tr ((£%£) (g, l n |V)M) , (32) 

where 

VjVtyVo = V^VSV Q V = l n (33) 

for any unitary Vo and V. Since the right-hand sides of (J3T|) and are equal, their 
left-hand sides must be equal as well, which concludes the proof of (J30J). Additionally, 
if an operation 17% weakly commutes with Ui, i.e., UiU^ ~ U^Ui, then L 7 ^ is an element 
of U(2) ®U(n). 

3.6. Comments and remarks 

Due to the arbitrary nature of the encoded component V (V G U(n)), every logical 
operation is represented in the (2n)-dimensional space by an infinite number of logically 
equivalent transformations = £ ® V (which collectively form an equivalence class 
characterized by £). Therefore, the use of MLE of logical basis states allows for a great 
flexibility in the choice of the actual physical transformation that realizes the desired 
logical operation, which can be an important advantage in many practical situations. 

An interesting question is whether the flexibility in the encoded component of 
makes it easier to find an optimal control field that produces a physical transformation 
belonging to the desired equivalence class of logical operations. The answer is not yet 
clear since £, the logical component of 17%, must still be produced with the same accuracy 
as in the single-level representation. Another question for future research is whether the 
flexibility due to the logical equivalence of different physical transformations improves 
the robustness of logical operations against noise in controls. 

A restriction in the MLE formalism presented above is the requirement that both 
logical basis states, |0) and |1), are encoded by the same number of physical levels. 
Satisfying this requirement in a real physical system can sometimes be a difficult task. 
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Consider for example, a logical basis encoded into vibrational levels on two electronic 
surfaces of a molecule. If the difference between the two vibrational frequencies is 
sufficiently large, a laser pulse of a given spectral width will encompass a different 
number of vibrational levels on the two surfaces. Therefore, it is necessary to consider 
how restricting the encoding subspaces to be of equal dimension can be dealt with 
in practice. First note that in any practical situation it is impossible to achieve the 
target quantum operation exactly. Correspondingly, laboratory controls are designed 
to produce an operation which is as close as possible to the target (an optimal control 
problem). In the case of SLE, the actual quantum operation should be as close as 
possible to the 2x2 target operation £. Analogously, in the case of MLE, the actual 
quantum operation should be as close as possible to the (2n) x (2n) target operation 
U$ = £ ® V (with the additional flexibility provided by the arbitrary nature of V). 
Consider now what will happen if in reality one logical basis state is encoded into n 
levels and the other into m = n + k levels {k > 0). In such a case the actual operation 
will be represented by an (n+m) x (n+m) unitary matrix C/i a b, and the target operation 
should be of the form 



where W is an arbitrary k x k unitary matrix (this mathematical structure also frequently 
appears in the context of QEC and operator QEC, e.g., see [I3J El)- The physical 
meaning of this target operation is that a qubit encoded in the (2n) -dimensional space 
should not couple to the additional k levels. So the laboratory task is to design controls 
which will produce the actual operation £/i a b as close as possible to target of (}34"j) . This is 
an optimal control problem which is not different in principle from the one encountered 
when n = m. Of course, the practical difficulty of finding optimal control fields can 
increase when the encoding subspaces have different dimensions, but the method of 
closed-loop laboratory control with shaped laser pulses is well-suited for dealing with 
such problems. Moreover, one can argue that optimal control could benefit from freedom 
in the choice of the actual parceling of the vibrational levels between the two subspaces 
since the experimenter does not have to a priori define the levels into pre-assigned 
groups. The most effective partitioning of physical levels into the logical basis can be 
included as part of the optimal control problem, with this additional freedom providing 
a potential for better solutions. 

4. Multilevel Encoding for Mixed States 

The principle of logical equivalence can be extended to mixed states. Consider a density 
matrix that represents a state in the expanded (2n)-dimensional Hilbert space: 



where are matrix elements representing the state of the logical component and Rij are 
n x n matrices representing the state of the encoded component. As a density matrix, 



target 



(f ® V) ® w, 



(34) 




(35) 
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p has certain fundamental properties: p = p\ (ip\p\ip) > for all |^>), and Tr (p) = 1. 
In Appendix A we show that the logical-component matrix 



Pl = ( ni n2 (36) 
r 2 i r 22 J 

as well as the diagonal sub-blocks Rn and R22 also have these properties and thus are 
proper density matrices. 

Logical equivalence for two density matrices p and p' is defined as 

Tr(pA r ) = Tr(p'A r ), r = {x,y,z}, (37) 

plus the normalization condition Tr (p) = Tr (p') = 1. Two unitary operations Ui and 
U2 are logically equivalent if they satisfy 

U lP Ul ~ U 2 pU\. (38) 

By evaluating the corresponding traces, it is straightforward to show that for the weak 
identity Uj, 

p ~ VipU\ (39) 

and 

UiU iP U\u\ ~ u^u lP u\u\ 

for any E U(2) ® U(n), (40) 

which is completely analogous to the result obtained for pure states in the previous 
section. 

Now we investigate the process of non-unitary evolution (e.g., due to 
environmentally-induced decoherence) in the context of MLE. We restrict our 
consideration to processes which are unitary in the logical component and general (not 
necessarily unitary) in the encoded component. Quantum operations (acting on the 
reduced density matrix p of the system of interest) representing such processes will 
be symbolically denoted as = {£,W}, where £, the logical component, is a 2 x 2 
unitary matrix and W, the encoded component, is a generalized quantum operation 
whose properties will be determined below. The action of £g on p is represented as 

2 

Using the explicit matrix form of this can be rewritten as 

2 

AP= E (fe»-i*6fcl»>Lu(i|)®aVi2,*). (42) 

i ,j,k,l=l 
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Since and W act on the density matrices p and {Ru}, respectively, these operations 
must preserve Hermiticity, positivity, and the trace. Therefore, a permissible operation 
W is a positive and trace-preserving map. 

Next, consider the effect of the generalized weak identity, Cj = {1 2 , W}, on p: 

2 

P = Cl P = E (^IOl lO'D ® • (43) 

ij=l 

Using the explicit form of the EA operators, 

2 

A r = A r <8> l n = ^ A^I^l L (j| ® l n , (44) 

we obtain 

2 

Tr(p'A r ) = ^r^A^TrCW^) 

2 

= ^ r^A^Tr (/?,,) = Tr(pA r ), (45) 

where we used the fact that W is a trace-preserving operation. This proves that Cj 
indeed has the property of the weak identity: 

Cip ~ p. (46) 

An immediate and interesting consequence of ()46|) is the logical equivalence of mixed and 
pure states. In other words, while p^ = \ip){ip\ represents a pure state in the form of (jlj), 
its logically equivalent state p' = Cip^ will be mixed if Cj includes any permissible and 
non-unitary operation W. Thus, a logical basis state can be encoded into a subspace 
of multiple physical levels which are mixed (e.g., by a thermal excitation) without any 
loss of coherence between the two logical basis states that form a qubit. For example, if 
an initial state is such that mixing is only within an encoding subspace and not between 
the two subspaces, then quantum computation with MLE can be conducted without 
any loss of fidelity due to intra-subspace mixing. This property makes the method of 
MLE a very promising practical alternative to the cumbersome process of ground-state 
cooling. 

An additional property of W is found by examining the weak commutation relation. 
Applying the method of proof used in (|31|) and (|32|) to general quantum operations in 
the form of (|41)1. it follows that any permissible operation weakly commutes with 
£j, i.e., 

dC^p ~ Ct&p, (47) 

if 

Tr (WoWRij) = Tr (WW Rij) for all i, j. (48) 
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Here, Wo and W are the encoded components of Cj and respectively. Equations (|T7|) 
and ()48|) show that all permissible operations {W} form a closed set. In other words, 
if W and W are two arbitrary permissible operations, then W" = WW is another 
permissible operation. An example of a positive and trace-preserving map with this 
closure property is the Kraus representation: 

WRij = W v RijW$, where ^ W*W V = 1„, (49) 

V V 

which is in fact a completely positive map 0. 

It is now possible to generalize the notion of a class of logically equivalent operations 
to non- unitary processes. For each logical component £ £ U(2), there exists an infinite 
number of equivalent logical operations (forming an equivalence class) which permit 
decoherence within an encoding subspace but preserve coherence between the encoding 
subspaces. Specifically, any two operations £g = {£, W} and £ = {£, W'}, with 
permissible maps W and W", are logically equivalent (also illustrating that a non-unitary 
operation £g = {£, W} is logically equivalent to a unitary operation C/g = £ £g> I 7 , where 
V £ U(n)). Accordingly, MLE allows for quantum computation with full fidelity in the 
presence of non-unitary evolution within the encoding subspaces. 

5. Example: Numerical simulation of closed-loop optimal control of unitary 
operations in a multilevel system 

To illustrate the formalism and advantages of MLE, as developed in sections 120 we 
perform a simulation of closed-loop optimal control of single-qubit unitary operations 
for a model system. It should be emphasized that the method of closed-loop optimal 
control is uniquely suited for laboratory implementation, but is very difficult to simulate 
on the computer, except for the simplest models. This is related to the fact that learning 
algorithms (e.g., genetic algorithms) which are used to search for the optimal control 
fields, involve vast amounts of data processing. In our numerical simulations, solving the 
time-dependent Schrodinger equation for each set of control parameters is the limiting 
computational step. Since a learning algorithm typically needs to search over hundreds 
of thousands or even millions of points in the parameter space until a global optimum is 
found, the simulation can take days even for relatively simple model systems. However, 
in the laboratory, a real physical system, whatever complex it may be, "solves" its own 
time-dependent Schrodinger equation in real time. Hence, the data processing of the 
learning algorithm in the laboratory is limited mainly by the repetition rate of the laser 
and pulse-shaping system, which can be as fast as several kilohertz and consequently, the 
search over a million of different control fields can be completed in just several minutes. 

5.1. The model 

Due to the previous considerations, we perform just a proof-of-principle numerical 
simulation of optimal control for a simple four-level model system which nevertheless 
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demonstrates the difference between the MLE and SLE operations in the presence of 
decoherence. This model represents a single qubit with each logical basis state encoded 
either by two levels in the case of MLE (n = 2) or by one level in the case of SLE (the 
other two levels remain unused). The energy spacings correspond to the electronic and 
vibrational levels of the sodium dimer, Na2-§ For MLE, the ground and excited electronic 
states correspond to the |0) and |1) logical basis states, respectively, i.e., each logical 
basis state is encoded by a pair of vibrational levels on the corresponding electronic 
surface. For SLE, the logical basis states |0) and |1) are defined as the lower vibrational 
levels on the ground and excited electronic surfaces, respectively. In a homonuclear 
dimer such as Na2, transitions between vibrational levels within a given electronic 
surface are forbidden (the dipole moment for these transitions is zero). Transitions 
between different electronic surfaces are allowed and correspond to the dipole moment 
set to unity. Thus, the number of allowed transitions equals n 2 . 

Unary multilevel atomic and molecular systems, such as Na2, are not scalable 
for quantum computing applications in the absence of physical entanglement. 
Implementations of quantum algorithms without entanglement have been investigated 
using Rydberg atoms [HJ EZj and linear optical systems |1H1 ED] • It was concluded 
that quantum computation with a single multilevel system is possible, but requires 
exponentially greater overhead than a multiparticle quantum computer that operates 
with entanglement |5T| E2] ■ Therefore, there exists significant interest in engineering 
interactions between atoms or molecules for creating entanglement needed for efficient 
and scalable quantum information processing [HH1 EH E3 EHJ EH EBJ EH E01 EH E2j- 
However, in our example scalability is not the issue, as we use the model multilevel 
system to investigate some of the most basic features of MLE for a single qubit. 

There exist plans to use homonuclear dimers such as Na2 and K 2 for experimentally 
studying the basic functioning of MLE in single-qubit systems in the presence of 
thermalization. These experiments would benefit from existing technologies for 
managing and measuring the dynamics of vibrational molecular wavepackets [23 ESI, 
including applications of closed- loop control methods [HI] . A much greater challenge 
would be an experimental realization of MLE in entangled multi-qubit systems, such as 
photon-atom quantum memories, as the technology for this future type of experiments 
is not yet fully developed. 



5.2. Formulation of the optimal control problem 

The time evolution operator U(t) for an isolated quantum system satisfies the 
Schrodinger equation: 

jU(t) = ~H(t)U(t). (50) 

§ The energy separation between the ground and first excited electronic surfaces is 0.066889653 a.u.; 
the separations between vibrational levels are 0.0007250238 a.u. and 0.000534563 a.u. on the ground 
and excited electronic surfaces, respectively. The data are taken from [35] , 



Encoding a qubit into multilevel subspaces 



17 



The time- dependent Hamiltonian, H(t), in this model is 

H(t) = H Q -fj£(t), (51) 

where Ho is the system Hamiltonian in the absence of control, /i is the electric dipole- 
moment operator and e(t) is the time-dependent control field, defined as 

n 2 

<t) = fit) a * cos ^ + 5 *)- ( 52 ) 

i=l 

Here, f(t) is an envelope function which incorporates the laser pulse width (e.g., a 
Gaussian or similar type of distribution) and at, u>i, 5i are the amplitude, frequency, and 
relative phase of the zth electric field component, respectively. Transition frequencies 
are determined by the system Hamiltonian. The operation is over the time interval 

[o,t/MI 

The control objective is to achieve the time evolution operator U(tf) which is as 
close as possible to the target transformation [/target- The fitness of the control field is 
evaluated by using the gate fidelity, which is a functional of control: 

^[e(t)] = lHl^ar g et-t/(t/)||. (53) 

Optimal control solutions correspond to maxima of T . The search for a global maximum 
is performed using a genetic algorithm (GA) implemented with population sizes of 
~ 200, several different reproductive schemes, and crossover and mutation rates between 
20 and 40 percent. 

5.2.1. Fidelity Junctionals for MLE. For MLE, the specific matrix norm in (J53j) will 
depend on the target transformation. The general form of the unitary evolution operator 
U(tf) can be written as 

U(t f ) '(el), (54) 

where A, B, C, and D are square matrices of dimension n. For the bit-flip operation U x 
of (f2T)J) and the phase-flip operation U z of (|2*3|) as the target transformations, we use 
fidelity functionals of the form 

(\\BB^-l\\ 2 ||£-C|| 2 \ 1/2 
^=l-(L^ + L_^j (55, 

and 

/||A4t-l|| 2 \\A + Df\" 2 , v 

* = i-(L_J- + J!_+JLj , (56) 

respectively, where the matrix norm is defined as ||M|| 2 = Tr (MM^ and the coefficients 
are chosen for the proper normalization, so that < T x < 1 and < T z < 1. The 
maximum of the fidelity functional [T x = 1 or T z = 1) is achieved when the actual 
transformation U(tf) is an element of the corresponding equivalence class (U x or U z ). 
|| In our simulations, tf can vary from 2 x 10 3 to 5 x 10 4 a.u. (i.e., from about 48 fs to 1.2 ps). 
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5.2.2. Fidelity junctionals for SLE. For SLE, two vibrational levels are used to encode 
the logical basis states, while the other two levels are unused. Therefore, the bit-flip 
and the phase-flip operations are represented by the unitary transformations 

U' x = a x @v (57) 

and 

U' z = a z ® w, (58) 

respectively, where a x and o~ z are the corresponding Pauli matrices acting on the space 
of the logical basis, and v and w are arbitrary 2x2 unitary matrices acting on the space 
of the two unused levels. For computations, we rearrange U' x and U' z so that their indices 
correspond to the physical order of the energy levels (the logical basis states |0) and |1) 
correspond to the 1st and 3rd energy levels, respectively). For the SLE operations U' x 
and U' z as the target transformations, the fidelity functional will be of the form 

^=5^13 + ^311 (59) 

and 

F M = -\Uu-U a \, (60) 

respectively, where Uy denotes the corresponding matrix element of the actual 
transformation U(tf). The maximum of the fidelity functional [T' x = 1 or T' z = 1) 
is achieved when the actual transformation U(tf) produces the target transformation 
(a x or <t 2 , up to a global phase) in the two- level subspace of the logical basis. 

5. 3. Numerical optimization and analysis of results 

Table ^ summarizes the gate fidelities obtained from the GA optimization for MLE 
(n = 2) and SLE configurations of the four-level model system. The GA is capable 
of finding optimal control solutions for both MLE and SLE cases with reasonable 
accuracies, limited mainly by the parameterized structure of the control fields in ([52)1 . 
These fidelities could be improved by coupling the results of the GA to an additional 
gradient-based search algorithm that does not rely on any particular parameterization 
of the control fields. However, for the purposes of our simulations, the logical operations 
obtained with the GA were sufficient to demonstrate the main distinctions between MLE 
and SLE. 

We are also interested to see what happens when we apply the optimal 
transformation located by GA optimization on an initial state which is mixed by 
a decoherence process (we assume that the effect of decoherence during the logical 
operation is negligible and that only the initial state is affected). We consider two types 
of decoherence processes: dephasing (only for MLE) and thermalization (for both MLE 
and SLE). 
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Table 1. Fidelities of MLE and SLE logical operations using the GA optimization. 







T' 


T' 


0.9987 


0.9975 


0.9998 


0.9996 



To investigate the effect of initial state decoherence on the fidelity of quantum 
operations, we consider the difference between the perfect (or intended) initial state 
rotated by the perfect target transformation: 

Ptargct = ^target pfAarget; (61) 

and the actual initial state rotated by the optimal actual transformation found by the 
GA: 

P(tf) = U(t f )pU(t f )K (62) 

The actual initial state p may be affected by a decoherence process and therefore differ 
from the perfect initial state p. For MLE, the difference between p t ar g et and p(i/) is 
estimated by evaluating the "error of equivalence" : 

£ iP> P\ = \ l Tr ( Ar M*/) ~ target]) I > (63) 

r=x,y,z 

where {A x , A y , A z } are the three EA operators of (J7J). This error measures how far the 
two density matrices are from complete logical equivalence. The worst value of e is 1, 
while in the case of complete equivalence e = 0. For SLE, the distance between the 
perfect target state and the actual final state is measured by 

e'[p,p] = Tr([p(t / )- Ptargct ] 2 ). (64) 

In the absence of decoherence (p = p), the errors e and e' will be non-zero only due to 
the fact that the optimal control solution found by the GA produces a transformation 
that is slightly different from the target one. However, if decoherence impairs the initial 
state (p 7^ p), this can significantly affect the error. We will see below that operations 
resulting from MLE are not at all hindered by dephasing or thermalization of the initial 
state. In contrast, the performance of SLE operations deteriorates if the initial state is 
thermalized. 

5.3.1. MLE operations in the presence of dephasing. Dephasing of non-stationary 
vibrational wave packets in hot alkaline dimers is caused by the vibrational-rotational 
coupling [T^]. The typical dephasing time is inversely proportional to the temperature 
and varies from ~ 3 ps for very hot molecules (produced in a heat pipe at about 400 
°C) to ~ 20 ps for vapours cooled to about 100 K. Note that the decoherence time 
on the electronic surfaces (i.e., within the encoding subspaces of the MLE scheme) due 
to this coupling is much shorter than that between electronic surfaces associated with 
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Table 2. Average errors of MLE operations with and without dephasing of the initial 
state. 



Initial state 


(e) for U x 


(e) for U z 


P 


0.0093 


0.0146 


P 


0.0076 


0.0134 



spontaneous emission from the excited electronic surface, which happens on the time- 
scale of several nanoseconds. Therefore, dephasing induced by the vibrational-rotational 
coupling is the most important mechanism of decoherence for vibrational wave-packets. 
Since the length of the optimal control pulse does not exceed 1.2 ps (and can even be 
made shorter at the cost of a small decrease in the gate fidelity), the effect of dephasing 
during the logical operation is negligible at temperatures about 100 K. 

We study the effect of dephasing on MLE operations by considering a set of 
randomly generated pure initial states of the form 



cosf^o) \ M , / cos(#i) \ 
exp (z0 o ) sin (0 O ) J y exp {tfa) sin (9 1 ) J 

and 9i are randomlv generated real parameters such 1 



where c , Ci, (fi , <pi, 9 , and 9\ are randomly generated real parameters such that 

Cq + c\ = 1, (66a) 
< 0o, 0i, 0o, 0i < 2vr. (66b) 

The initial state l^) corresponds to the density matrix p = \if))(if)\, whose elements in the 
basis of the four energy levels are denoted as (i\p\j) = = 1,2,3,4). Dephasing 

within each encoding subspace is represented by setting the off-diagonal matrix elements 
of each diagonal sub-block equal to zero, resulting in the mixed state 



P 



( r n r 13 r M \ 

r 22 r 23 r 24 

r 3 i r 32 r 33 

y r 4 i r 42 r 44 J 



(67) 



The error of equivalence for the U x and U z target transformations is evaluated both in 
the absence and in the presence of dephasing: e[p, p] and e[p, p], as defined in (jHHjl. and 
averaged over 10 6 random choices of the initial state in the form of (|65j) . These 
average errors, (e), are presented in table El We see that the operation error does not 
increase (and even slightly decreases) when the initial state is affected by dephasing. 
This means that the pure state p and the dephased state p are equally suitable for MLE 
operations due to their logical equivalence. 



5.3.2. MLE and SLE operations in the presence of thermalization. As mentioned 
above, alkaline dimers are typically produced at high temperatures of up to 400 °C. 
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Moreover, these molecules are unstable at low temperatures, although vapours of 
sufficiently low concentrations can be carefully cooled to 100 K or even slightly below. 
We study the effect of thermalization on both MLE and SLE operations by considering 
the following pure initial state: 

m = |0)l ® ( J j =(1,0,0, Of, (68) 

which corresponds to the lowest energy state of the molecule. The resulting density 
matrix p = I^K^I nas onr y one non-zero element: (i\p\j) = SuSji = 1,2,3,4). 
Thermalization impairs the pure initial state of (|HH)l and results in a mixed state p of 
the form 



/ 1 


-A 








o\ 







A 






















V 














(69) 



where A = exp(—E v /kBT), ks is the Boltzmann constant, T is the temperature, and E v 
is the energy separation between the vibrational levels on the ground electronic surface. 
Thermal excitations to all energy levels higher than \n v = 1) are neglected, which is a 
reasonable approximation at temperatures of about 100 K and below. 

At zero temperature (i.e., when p = p), the errors of MLE and SLE bit-flip 
operations (as defined in (|6l?|) and (ffiljl. respectively) are e = 0.0097 (for U x ) and 
e' = 0.0009 (for U' x ), respectively. This is consistent with the fact that optimal controls 
found by the GA have a higher fidelity for the SLE operations in comparison with the 
MLE ones, as shown in table Q However, the situation is strikingly different in the 
presence of thermalization. We evaluate the errors for both MLE and SLE operations, 
e[p, p] and e'[p,p], respectively, for the range of temperatures between 70 K and 120 K 
(corresponding to values of A between approximately 0.037 and 0.129). These errors 
are shown in figure El versus the temperature. The error of the SLE operation increases 
quite rapidly with the temperature and at temperatures above 85 K the performance 
of the SLE operation is worse than that of the MLE one. This increase of the SLE 
operation error is explained by the fact that at higher temperatures a larger portion 
of the population is transferred out of the logical basis. On the other hand, the 
performance of the MLE operation changes very little with thermalization (in fact, the 
error of the MLE operations slightly decreases as the temperature increases). This is 
explained by the fact that the population transfer within the encoding subspace (caused 
by thermalization) does not affect the logical state of a qubit with MLE, as discussed 
in section HI Thus, the pure state p and the thermalized state p are equally suitable for 
the MLE operation due to their logical equivalence. 
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Figure 2. Errors of MLE and SLE bit-flip operations, e[p,p] of (jfiHfl and e'[p,p] of 
respectively, for the thermalized initial state p of (ff!§)l . versus the temperature. 



6. Two-Qubit Operations for Multilevel Encoding 

For a two-qubit system, the formalism of MLE is analogous to that developed in 
the previous sections for a single qubit. As before, the energy levels of a quantum 
system realizing the kth qubit are divided into two distinct encoding subspaces of equal 
dimension, Sok and Sik- Let \a)k denote any element of Sok such that k(a,\a)k = 1 and 
\b)k denote any element of Sik such that = 1, i-e., 

T 



I a )k — a ik\Xi)k — [o-lk a nkj 

i=l 

n 

such that \aik\ 2 = 1 (70) 

i=l 

and 

n rp 

\b)k = ^2b ik \(pi)k = (bik ■ ■ ■ b nk j 

n 

such that \bik\ 2 = 1, (71) 

i=l 

where k G {1, 2} denotes qubit 1 or qubit 2, n is the dimension of all encoding subspaces, 
and and {\4>i)k} are orthonormal bases that span Sok and <Sifc, respectively. 
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Correspondingly, k {Xi\Xj)k' = k{H4>j)k' = <M**' and k{Xi\<t>j)k< = 0. Now the logical 
basis of the kth. qubit is expressed as 

|0) fc = |0) L ® |o) fc = (a lk ■■■ a nk (72a) 

n 

\i) k = |i) L ® |&)* = (^_o hk ■■■ Kk) T - (72b) 

n 

It follows that fe(0|0)fc' = fc(l|l)fe' = S kk > and fc(0|l)fc/ = 0. In the abbreviated notation, 

io> - ("*»')• ii)i= U)- (73) 

and the state vector for the kth qubit can be written as 

|^ = c ofc |0> fe + c lfc |l> fe = ^* ) > (74) 

where |cofc| 2 + |cifc| 2 = 1 (to ensure normalization) and a k = co k \a) k and (5 k = c\ k \b) k . 
The logical basis of the two-qubit system can be defined as 

|0) 1 ®|0> 2 =(' |a >1 U(' |a )2 Y (75a) 

l )^®! 1 )^ ( '"o 1 ) ® ( |6°> 2 ) ' (75b) 
W^W^f IJ), J ® ( ' a >2 j. (75c) 

|1>lS>|1> ^(l^) S> (l^)- (7M) 

However, it is much more convenient to use another basis, which is obtained by a unitary 
transformation on the basis defined above. The preferred logical basis for the two-qubit 
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(TQ) system is 



|l>TQ = i2(|0>l® |0); 



|2>tq = 72(10)! ®|1) 2 ) 



( \a)x® | o) 2 \ 




\ 

|a)i <g> 1 6)2 





(76a) 



(76b) 



|3) 



TQ 



i2(|l>i® |0>! 



|4) TQ = i2(|l)i<8) |1) 2 ) 






|6)i ® \a) 2 


\ 



V 16)1 ® 1 6)2 y 



(76c) 



(76d) 



where is a 4n 2 x 4n 2 unitary transformation. Now consider a state vector of the 
two-qubit system. In the logical basis of (fTTjj) . \^/) is expressed as 



i*> = X>io 



TQ, 



(77) 



i=l 



where 5^i=i |cj| 2 = 1 ensures that (^l^) = 1 (|^) is a vector of dimension 4n 2 ). 

In order to formulate the principle of logical equivalence for the two-qubit system, 
we introduce a set of 15 EA operators, 



A r — A r 



1, 



r = {l,2,...,15}, 



(7f 



where l n 2 is the n 2 x n 2 identity matrix and {A r } is the set of 4 x 4 matrices which 
represent the 15 generators of SU(4), which, with U 4 , span the space of all 4 x 4 matrices. 
The principle of logical equivalence is formulated analogously to the single-qubit case. 
Two physical states, |^) and are logically equivalent (denoted as |^) ~ |^')) ^ 



(*|A r |*> = (* , |A r |* / ), r = {l,2,...,15}. 



(79) 



Now we develop the explicit form of logical operations for the two-qubit system with 
MLE of the logical basis states. First, consider local operations which act separately on 
the two qubits. Such an operation is represented (in the preferred basis) by the matrix 



U = R ([/« <g> U^) E}, 



(80) 
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where acts on the kth qubit. In particular, the weak identity operation that 

maintains logical equivalence for the entire two-qubit system and simultaneously for 
each of the individual qubits is given by 



Ui = R(ujp ®Ul 2) )R* = 1 



(2) 



(1) 



Vr 



(2) 



where C/j is the weak identity operation and Vq is an arbitrary n x n unitary matrix 
acting on the kth qubit. As in the single-qubit case, the two-qubit weak identity satisfies 

17/1*} ~ |*) for any (82) 

The proof given in (}3Tj) and (|52^1 easily extends to the two-qubit case, so any 
operation U such that 



U e U(4) <g> [U(n) <g> U(n)] C U(4n 2 



(83) 



satisfies UjU ~ UUj. Therefore, any two-qubit logical operation S resulting from SLE 
will be represented in the framework of MLE by a matrix of the form 



Us 



S <g> [V® v^] 



(84) 



where S G U(4) and e U(n). Since and are arbitrary unitary matrices, 
any logical operation will be represented by an infinite number of logically equivalent 
physical transformations. 

The tensor-product form ® of ()84j) is required in order to ensure that 
entanglement created by a logical operation is reversible not only by a specific inverse 
operation, but by any member of a class of logically equivalent operations. To 
demonstrate this point, consider the C-NOT gate, Uqn- For MLE, this operation is 
represented by the following set of logically equivalent {An 2 ) x (4n 2 ) matrices: 



U, 



CN 



( 1 \ 

10 

1 

x 1 J 

( ® V® 

® V® 



\ o o 
















55) 



/ 



Equation (J85|) defines an entire class of logically equivalent operations. With SLE, 
applying the C-NOT operation twice results in the identity operation. Thus, applying 
the operation {Tcn twice should result in a logically equivalent state, i.e., 



(86) 
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Table 3. The group-theoretic structure of unitary logical operations for single- level 
and multilevel encoding of qubits. 



Encoding 


Single qubit 


Two qubits 


M qubits 


n = 1 


U(2) 


U(4) 


U(2 M ) 


n > 1 


U(2) ® U(n) 


U(4) <g> [U(n) <8> U(n)] U(2 M ) 


8> [U(n) ® • • • <g> U(n)] 








M times 



Indeed it is clear that Uq n ~ Ui. Moreover, applying C/qn twice actually means that 
we can sequentially apply two physically different transformations which belong to the 
same equivalence class defined by (J85|) since the TPS assures that for any pair of logically 
equivalent operations U\ G Uqn and U 2 G Uqn, the overall operation U-JJi is equivalent 
to the weak identity Uj of (|ST|). 

The C-NOT gate of (|85|) . together with the Hadamard gate of (|26|) and phase-shift 
gate of (|27p. establish universality [S3] for quantum computation with MLE. 

7. Multi-qubit Operations for Multilevel Encoding 

In the case of M qubits, logical equivalence is based on the set of Z M = (2 2M — 1) EA 
operators 

A r = A r ®l n M, r = {1,2,..., Z M } (87) 

where l n M is the n M x n M identity matrix and {A r } is the set of 2 M x 2 M matrices which 
represent the Zm generators of SU(2 M ). A possible choice of these generators is, for 
example, the standard Cartan- Weyl basis or its Hermitian variant [66J. The principle 
of logical equivalence for a pair of multi-qubit states and |^') will be expressed as 
usual: 

~ |^') <<=> (^|A r |*) = (^'|A r |^') for all r. (88) 

In addition, the states must be normalized. 

Unitary logical operations for M qubits, with n-dimensional encoding of every 
logical basis state, will be elements of the group 

U(2 M ) ® [U(n) ® • • • ® U(n)l C U ((2n) M ) . (89) 

V v ' 

M times 

Table El summarizes the group-theoretic structure of this mapping from single-level to 
multilevel encoding. 

The generalization of the MLE formalism to mixed states and non-unitary logical 
operations for multiple qubits is straightforward and analogous to the single-qubit case 
(see section 0J . 
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8. Conclusions 

The formalism of MLE of logical basis states presented in this paper is motivated 
primarily by practical realities such as imperfections and noise in controls, thermally 
excited initial states, and environment ally- induced decoherence. 

With MLE, a given logical operation is realized by an infinite number of 
equivalent physical transformations which may be used interchangeably throughout the 
computation, thereby allowing for a significant flexibility in the laboratory realizations 
of control. A crucial point is that MLE naturally suits the application of ultrafast 
broadband controls which simultaneously drive multiple transitions and have the 
advantage of generating faster quantum operations, helping to lessen the effect of 
decoherence. 

As presented in section 03 with a simplified model of a single qubit based on Na2 
(which is not efficiently scalable for quantum computing applications in the absence 
of entanglement), the simulation of closed-loop optimal control of single-qubit MLE 
operations demonstrates that ultrafast optimal solutions for multilevel systems are 
readily achievable. Moreover, these solutions are not affected by various decohering 
processes which are detrimental to quantum information processing based on SLE. 

An open question is whether the multiplicity of logically equivalent physical 
transformations will facilitate the discovery of controls with increased robustness against 
noise in the control fields. While this assumption is intuitively plausible, a rigorous 
analysis of the robustness is still required to make any formal conclusions. A powerful 
method for the robustness analysis is the study of the landscape for optimal control 
solutions [HZ1 EH EH] ■ Other numerical approaches to the robustness analysis are being 
considered as well. 

An extremely promising feature of MLE is the ability to work with mixed initial 
states and decoherence without a loss of operational fidelity. If mixing and decoherence 
are contained within the encoding subspaces, then coherence between the logical basis 
states will not be affected, meaning that quantum computation with full fidelity is 
possible. This property, illustrated by a numerical example in section makes MLE a 
very attractive approach to practical quantum computation. 

The tensor-product structure of the Hilbert space, which appears in MLE, is a 
very general feature of multiparticle systems. Therefore, there are many possibilities for 
building an MLE structure in various physical systems. Interesting examples of systems 
which could be suitable for MLE are molecules and trapped ions (with encoding of the 
logical basis into subspaces of vibrational levels of a molecule or motional sidebands of an 
ion, respectively). Although we are planning experiments involving diatomic molecules 
such as Na2 and K 2 to study the basic functioning of MLE in ultrafast single-qubit 
operations in the presence of thermalization, we cannot yet specify particular quantum 
information systems which will benefit most from the use of MLE, since this depends 
on many technical considerations. Nevertheless, we estimate that the tendency toward 
faster quantum operations will definitely favour the application of MLE. 
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Appendix A. Properties of the Density Matrix Components 

In (135)1 . we expressed the density matrix of an arbitrary state in the expanded (2n) 
dimensional Hilbert space as 



<.=E'«i'Vki®*« = ( r °* a ), (a.i; 

~~ x \ T2\IV2.\ r 22-K22 ' 



where are matrix elements representing the state of the logical component and Rij are 
n x n matrices representing the state of the encoded component. This is a generalization 
of a pure state, (jU): 

|V>)=co|0)l®|o) + ci|1) l ®|6). (A.2) 

If the dimension of the encoding subspaces decreases from n to 1, then p and \ip) will 
represent the logical state of a two-dimensional quantum system: 

Pl = ( r " ) (A.3) 
\ r 21 r 22 J 

and 

|^)l = c |0) l + Ci |1) l . (A.4) 

Thus, for encoding subspaces of any dimension, the logical component (pl or 
represents the state of a two-level qubit. Therefore, pl is a proper density matrix, i.e, 
it is Hermitian, positive, and normalized. This normalization means that Tr (pl) = 
+ r 22 = 1, where r n , r 22 > 0. 
In this appendix, we will also show that the matrices Rn and R 22 of the encoded 
component are proper density matrices that satisfy all necessary properties. The 
Hermiticity of Rn and R 22 follows directly from the Hermiticity of p. Examining the 
matrix structure of p in (jA.ljl . consider the following scenarios: 

(i) If 1 < k, I < n, then p k[ = (Rn)ki, which implies that (Rn)ki = (Rn)iki meaning 
Rn = R\i- 

(ii) If n + 1 < A;, I < 2n, then p ki = {R 22 ) k 'v, which implies that (R 22 ) k 'v = (-^22)*/^, 
meaning R 22 = R\ 2 (where k' = k — n and V = I — n). 

However, if 1 < k < n and n + 1 < / < 2n, then p k i = (Ri 2 ) k i>, which implies that 
{Ri2)ki> = (#21) iw, meaning R 12 = R) 21 . 

The positivity of Rn and R 22 also follows directly from the positivity of p: 

(V>l/#> > for all (A.5) 
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where 

= \c \ 2 r n (a\R n \a) + clc^uialR^lb) 

+ c c[r 2l {b\R 21 \a) + | 1 2 r* 22 < & | J^ 22 1 a> > 0. (A.6) 

In particular, if \ip) = |0)l <S> \a), then (ij;\p\ij;) = r 11 (a|i? 11 |a) > 0, which implies that 
(a\Rn\a) > for any |o). Similarly, if = \b), then (b\R 22 \b) > for any |&). 

Now consider the normalization condition for the density matrix p: 

Tr (p) = r n Tr (R n ) + r 22 Tr (R 22 ) = 1. (A.7) 

Using the normalization of the logical component, rn + ri 2 = 1, we obtain 

r n Tr (i2 u ) + (1 - r n ) Tr (i? 22 ) = 1, (A.8) 

Since the logical and encoding components are independent, i.e., any logical state can 
be "attached" to any encoding configuration, (jA.8|) is satisfied only when 

Tr( J R 11 )=Tr(i? 22 ) = l. (A.9) 

However, note that Tr (R12) 7^ 1 and Tr (R21) 7^ 1 in general. 

In conclusion, this appendix demonstrates that the density matrix components p^, 
Rn, and i? 22 are Hermitian, positive, and normalized (of unit trace), and therefore are 
proper density matrices. 
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